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ABSTRACT 

We report on a sensitive survey for radio pulsar wind nebulae (PWN) towards 
27 energetic and/or high velocity pulsars. Observations were carried out at 1.4 GHz 
using the Very Large Array and the Australia Telescope Compact Array, and utilised 
pulsar-gating to search for off-pulse emission. These observing parameters resulted in a 
considerably more sensitive search than previous surveys, and could detect PWN over a 
much wider range of spatial scales (and hence ambient densities and pulsar velocities). 
However, no emission clearly corresponding to a PWN was discovered. Based on these 
non-detections we argue that the young and energetic pulsars in our sample have 
winds typical of young pulsars, but produce unobservable PWN because they reside 
in low density (n ~ 0.003 cm -3 ) regions of the ISM. However, non-detections of PWN 
around older and less energetic pulsars can only be explained if the radio luminosity 
of their winds is less than 10 of their spin-down luminosity, implying an efficiency 
at least an order of magnitude smaller than that seen for young pulsars. 
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1 INTRODUCTION 

Almost all radio pulsars have rotational periods which are 
steadily increasing with time. This spin down corresponds 
to a loss of rotational kinetic energy E = 4tt 2 IP I P 3 , where 
I is the moment of inertia of the neutron star (assumed to 
be 10 45 g cm 2 ) and P is its period; for the known pulsar 
population E falls in the range 10 28 — 10 39 erg s _1 . The 
bolometric luminosity of the radio pulses themselves is in 
all cases a tiny fraction of the spin-down luminosity E, and 
it is thought that most of the pulsar's spin-down luminosity 
is dissipated via a magneti sed wind populated by relativis- 
tic electrons and positrons (Rees & Gunn 1974; Michel 1982 



Kennel & Coroniti 1984). However it is still not well under- 



stood how this wind is produced, how it evolves as it flows 
away from the pulsar, what it is composed of, how its prop- 



* Hubble Fellow 

t Current address: Physics Department, Furman University, 3300 
Pointsett Hwy, Greenville, SC 29613, USA 



pulsars: general - radio continuum: ISM - stars: winds 



erties depend on those of the pulsar itself, or how it changes 
as the pulsar ages. 

Particles in the wind move along magnetic field lines 
as they stream away from the pulsar magnetosphere, and 
produce no observable emission. At some distance from the 
pulsar, the pressure of the wind is eventually balanced by 
an external pressure, and the resulting shock randomises 
the pitch angles of the relativistic particles. These particles 
consequently gyrate in the local magnetic field and produce 
synchrotron emission. The properties of the resulting pulsar 
wind nebula (PWN) can then be used to determine various 
parameters of the pulsar wind which produces it. At ra- 
dio wavelengths, PWN are characterised by an amorphous 
or filled-centre morphology, a moderate degree of linear po- 
larization (~20%), and relatively flat spectra (a ~ —0.3, 
S u oc v a ) (jWeiler & Panagia 1978j). 



Various types of PWN are produced, depending on the 
source of confinement of the wind. Young pulsars are of- 
ten still located inside their associated supernova remnants 
(SNRs), and the hot gas produced by the SNR blast- wave 
provides the confining pressure. These PWN, also known as 
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"plerions" , are typified by the Crab Nebula. If the SNR has 
dissipated, the confining pressure is then that of the ambi- 
ent interstellar medium (ISM). This results in much larger 
"ghost remnants", which hav e been proposed bu t not ob- 
served (Blandford et al. 1973; Cohen et al. 1983). In cases 



where a pulsar has a high space velocity, the ram pressure 
resulting from its motion can dominate the ambient gas pres- 
sure, resulting in a bow-shock PWN (e.g. Frail & Kulkarni 
1991). 

While PWN can tell us much about pulsar winds, the 
number of sources we have to study is small - at radio 
wavelengths, fewer than 10 pulsars have observable PWN. 
All of these pulsars are very young and have high values 
of E. While ~20% of SNRs have a "plerionic" component, 
in most cases the associated pulsar has not been detected, 
and without knowing the pulsar parameters it is difficult to 
constrain the properties of the PWN and the corresponding 
wind. Thus it is of considerable interest to target candidate 
pulsars with the intention of either finding new PWN, or 
determining upper limits on such emission. 

Many searches for radio PWN around energetic or fast- 
moving pulsars have been carried out, at varying resolutions 
and surface-brightness sensitivities, but usually with no suc- 
cess (e.g. Schonhardt 1974; Weiler, Goss & Schwarz 1974; 
Cohen et al. 1983). The most recent and comprehensive of 
these searches, and the only one to specifically target young, 
energetic or high velocity pulsa rs, was the recent survey of 
Frail & Scharringhausen (|l997|), hereafter FS97. FS97 im- 



aged regions around 35 pulsars with the Very Large Array 
(VLA) at 8.4 GHz, and found no nebular emission associ- 
ated with any of their targets. Their stringent upper limits 
allowed them to conclude that most pulsars put less than 
10" 6 of E into radio emission from a PWN, ~ 100 times 
less than observed for the young, high E, pulsars which 
power detected radio nebulae. Based on this result, FS97 
concluded that pulsar winds change in some way as pulsars 
age and slow down, such that they are no longer efficient at 
producing radio emission. 

However, despite these apparently constraining limits, 
in hindsight the observing parameters for this search were 
probably not ideal for looking for PWN. First, FS97 argued 
that for ambient densities of 1 cm -3 and pulsar velocities 
of 150 km s _1 , PWN around almost all their sources were 
likely to be unresolved, even at their high spatial resolu- 
tion of 0'.'8. However, other choices of ambient density and 
pulsar velocity can produce PWN with much larger angu- 
lar extents, resulting in a flux density limit much poorer 
than estimated by FS97. Furthermore, the maximum scale 
to which FS97 were sensitive was only 20"; they could not 
detect PWN larger than this at any flux density. Secondly, 
in 40% of their sample FS97 detected a point source at the 
position of the pulsar, but had no way of distinguishing be- 
tween the compact PWN they were looking for and the pul- 
sars themselves. By extrapolating pulsar flux densities from 
much lower frequencies, FS97 argued that the flux densi- 
ties they were detecting at 8.4 GHz were consistent with 
the expected pulsed fluxes, and hence concluded that they 
were not detecting any PWN in these data. Finally, a typi- 
cal PWN (spectral index a — —0.3) has a significantly lower 
flux density at their observing frequency of 8.4 GHz than at 
lower frequencies. 

Motivated by these points, we have undertaken an ex- 



tensive survey for PWN at both northern and southern de- 
clinations, with observing parameters chosen to give much 
greater sensitivity to PWN. First, we have observed at 
1.4 GHz, at which frequency PWN can be expected to be 
~70% brighter than at 8.4 GHz. Secondly, we have observed 
at a reduced spatial resolution of ~12", to give better sensi- 
tivity to extended structure. Thirdly, our observations are in 
telescope configurations whose shortest spacings correspond 
to a spatial scale of many arcmin. Finally, and most impor- 
tantly, all our observations have employed pulsar-gating, in 
which images are made from data taken only when the pulsar 
is off; we are thus sensitive to compact and unresolved PWN, 
which might otherwise be masked by the pulsars themselves. 

Initial results from this survey, using the Australia 
Telescope Compact Array (ATCA), have been presented 
in two previous papers. Gaensler et al. (1998), hereafter 
GSFJ98, reported the discovery of a faint PWN associated 
with P SR B0906-49, while Stappers, Gaensler & Johnston 
hereafter SGJ99, presented non-detections of PWN 



towards four pulsars. We here report on the remainder of this 
survey, consisting of 1.4 GHz pulsar-gated observations of 27 
more pulsars, using the VLA and ATCA. In we describe 
our observations and analysis, while in ^ we present non- 
detections of PWN towards these sources. In 94| we quantify 
the improvement in sensitivity of the current survey, and 
discuss the constraints we can put on the radio luminosities 
of pulsar winds from our data. 



2 OBSERVATIONS AND REDUCTION 

The 27 pulsars observed are listed in Table all were cho- 
sen for their high E and/or space velocit y. From this sample 
22 pulsars we re observed with the VLA (Napier, Thompson 
fc Ekers 198S ) , w hile the remaining 5 pulsars were observed 
with the ATCA ( [Frater, Brooks fc Whiteoak 1992| ). All ob- 
servations were carried out at frequencies near 1.4 GHz. 

The VLA observations were made in the C configura- 
tion, using a bandwidth of 25 MHz for 00 h < RA < 12 h and 
12.5 MHz otherwise, and with a phase centre correspond- 
ing to the catalogued pulsar position. The observing time 
for each pulsar was typically 15 min. Amplitudes were cal- 
ibrated using observations of 3C 286 and 3C 48, assuming 
1.4 GHz flux densities of 14.9 Jy and 16.3 Jy respectively 
(where 1 Jy = 10 -26 W m" 2 Hz" 1 ). ATCA observations 
were made in the 6C configuration, using a bandwidth of 
128 MHz (further subdivided into 32 spectral channels); am- 
plitude calibration was carried out using PKS B1934-638 
and assuming a 1.4 GHz flux density of 14.9 Jy. For ATCA 
observations, the phase centre was offset from the pulsar's 
catalogued position by ~ 1'; each pulsar was observed for 
approximately 12 hr. Antenna gains and instrumental po- 
larization were calibrated using observations of strong unre- 
solved sources in the vicinity of each pulsar; all four Stokes 
parameters were recorded. 

All observations were gated at the pulsar period in order 
to look for off-pulse emission at the pulsar position, using 
ephemerides supplied by A.G. Lyne. For the VLA, gating 
was carried out by phasing up the array on a nearby cali- 
brator, then integrating on the pulsar for a few minutes. The 
analogue sum of the signals from all antennas was formed 
from these data, and then folded at the apparent pulse pe- 
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Table 1. Pulsars surveyed for PWN. Uncertainties in the last digit are indicated (or are omitted when the uncertainty is less than one in 
the last digit); positional uncertainties do not include systematic errors due to calibration. Three pulsars were not detected, for reasons 
discussed in the text. PSR B1706— 44 was included as a test case. 



Pulsar 


Date 




Telescope 


Pulsar Position (J2000) 


Pulsar Flux Density 


Off-pulse RMS 


^min 


$max 




Observed 




RA 


Dec 


at 1.4 GHz (mjy) 


(mjy beam -1 ) 


(") 


(') 


B0114+58 


1999 


Jan 05 


VLA 








0.47 


19 x 10 


15 


B0136+57 


1999 


Jan 05 


VLA 


01:39:19.76(1) 


+58:14:31.68(7) 


4.8(4) 


0.30 


18 x 10 


15 


B0355+54 


1999 


Jan 05 


VLA 


03:58:53.680(5) 


+54:13:13.63(5) 


11(2) 


0.26 


15 x 10 


15 


J0538+2817 


1999 


Jan 05 


VLA 


05:38:25.090 


+28:17:09.41(1) 


4.1(2) 


0.40 


12 x 11 


15 


B0540+23 


1999 


Jan 05 


VLA 


05:43:09.665 


+23:29:06.167(1) 


32.1(1) 


0.42 


12 x 11 


15 


B0611+22 


1999 


Jan 05 


VLA 


06:14:17.020 


+22:29:56.680(5) 


6.9(4) 


0.40 


12 x 11 


15 


J0631+1036 


1999 


Jan 05 


VLA 








1.8 


13 x 12 


15 


B0656+14 


1999 


Jan 05 


VLA 


06:59:48.13(2) 


+14:14:21.0(4) 


1.5(2) 


0.26 


12 x 12 


15 


B0736-40 


1999 


Feb 21 


ATCA 


07:38:32.342 


-40:42:40.16 


87(5) 


0.19 


17 x 8 


4.5 


B0740-28 


1999 


Jan 05 


VLA 


07:42:49.038 


-28:22:43.331(2) 


25.7(2) 


0.53 


30 x 10 


15 


B1356-60 


1999 


Mar 10 


ATCA 


13:59:58.5(6) 


-60:38:07.73(9) 


12.5(3) 


0.29 


16 x 10 


4.5 


B1449-64 


1999 


Feb 19 


ATCA 


14:53:32.712(1) 


-64:13:15.51(1) 


19.5(5) 


0.07 


11x9 


4.5 


B1508-57 


1999 


Feb 18 


ATCA 


15:12:43.041(3) 


-57:59:59.94(3) 


7.3(7) 


0.21 


10 x 9 


1.8 


B1634-45 


1999 


Feb 21 


ATCA 


16:37:58.729(9) 


-45:53:26.7(2) 


1.4(3) 


0.16 


15 x 9 


4.5 


B1706-16 


1999 


Feb 02 


VLA 


17:09:26.44(1) 


-16:40:57.4(3) 


6.6(3) 


0.50 


20 x 12 


15 


B 1706-44 


1999 


Feb 02 


VLA 


17:09:42.52(2) 


-44:29:06(1) 


7.3(7) 


0.36 


91 x 10 


15 


B1718-35 


1999 


Feb 02 


VLA 


17:21:32.71(3) 


-35:32:46.4(7) 


5(5) 


1.6 


51 x 11 


15 


B1719-37 


1999 


Feb 02 


VLA 


17:22:59.04(7) 


-37:11:57(2) 


1.2(2) 


0.53 


51 x 11 


15 


B1727-33 


1999 


Feb 02 


VLA 








1.6 


40 x 9 


15 


B 1730-37 


1999 


Feb 02 


VLA 


17:33:26.74(3) 


-37:16:56(1) 


2(1) 


1.2 


45 x 12 


15 


B 1754-24 


1999 


Feb 02 


VLA 


17:57:29.371(1) 


-24:22:02.22(2) 


2.8(1) 


0.15 


25 x 10 


15 


B1821-19 


1999 


Feb 02 


VLA 


18:24:00.460 


-19:45:53.571(2) 


8.1(3) 


0.19 


21 x 11 


15 


B1823-13 


1999 


Feb 02 


VLA 


18:26:13.16(3) 


-13:34:49.9(7) 


2.2(2) 


0.32 


16 x 12 


15 


J1835-1106 


1999 


Feb 02 


VLA 


18:35:18.32(2) 


-11:06:16.6(4) 


3(2) 


0.58 


16 x 12 


15 


B1930+22 


1999 


Feb 02 


VLA 


19:32:22.63(1) 


+22:20:53.3(2) 


1.6(2) 


0.54 


12 x 11 


15 


B1933+16 


1999 


Feb 02 


VLA 


19:35:47.830 


+16:16:39.806(1) 


75(1) 


0.34 


13 x 11 


15 


B2011+38 


1999 


Feb 02 


VLA 


20:13:10.341(1) 


+38:45:43.225(7) 


11(1) 


0.82 


12 x 10 


15 


B2148+63 


1999 


Feb 02 


VLA 


21:49:58.71(3) 


+63:29:44.9(2) 


2.9(2) 


0.59 


14 x 10 


15 



riod to give an un-dedispersed pulse profile. A gate was then 
set on-line, such that one IF recorded on-pulse data while the 
other recorded off-pulse data, effectively giving two bins of 
possibly uneven size. The smearing due to dispersion across 
the band was sufficiently small for all pulsars that it was al- 
ways possible to completely separate on- and off-pulse emis- 
sion when choosing the gate. For the ATCA, visibilities were 
recorded at high time-resolution (typically 32 bins per pe- 
riod), and then folded at the apparent pulse period before 
being written to disk. Dedispersion (of 32 channels across 
the 128 MHz bandwidth) was carried out during data re- 
duction, and appropriate phase bins were then chosen to 
generate on- and off-pulse images. 

Data were edited and calibrated using the HIRIAD and 



MP S packages according to s tandard procedures (Greisen 
199 6;; flault fc Killeen 1998| ). On- and off-pulse images 



of a field containing each pulsar were formed using uni- 
form weighting. Each image was deconvolved using either 
the CLEAN algorithm (for fields containing primarily point 
sources) or a maximum entropy algorithm (for fields con- 
taining significant extended emission), and then smoothed 
with a Gaussian restoring beam. For some sources the region 
was significantly confused by extended structure; in these 
cases, we constrained the deconvolution process by gener- 
ating CLEAN boxes from lower resolution Molonglo Galactic 
Plane Survey (MGPS) or Northern VLA Sky Survey (NVSS) 



data ( |Green et al. 1999^ |Condon et al. 1998|) 



For each source, the position of the pulsar was deter- 
mined by fitting in the u — v plane to the difference of the 
on- and off-pulse data. For VLA data, a flux density for each 



pulsar was determined by measuring a flux in the on-pulse 
image, then scaling using the width of the gate used. For 
ATCA data, flux densities were measured directly from un- 
gated data. The sensitivity of each image was determined by 
measuring the off-pulse RMS at the pulsar position in each 
case. 



3 RESULTS 

The positions, fluxes and off-pulse sensitivities for the 27 
pulsars observed are listed in Table [l], along with the range 
of spatial scales to which each image was sensitive. In most 
cases, the pulsar was clearly detected in the on-pulse image, 
and was completely gated out in the off-pulse image. 

As a test of our sensitivity we included PSR B1706- 
44 in our VLA observations, a pulsar which is known to be 



embe dded in a candidate PWN (Frail, Goss & Whiteoak 



1994). This nebula was easily detected in our image, with a 



surface brightness and spatial extent similar to that obtained 
in previous data. 

All but three of the pulsars in our sample were detected 
in these observations. For detections the measured inter- 
ferometric position was compa red with the timing position 
given in the pulsar catalogue ( jTaylor, Manchester fc Lyne 



199J). Our gated position for PSR B 1706- 44 differs signif- 



icant ly from the original timing position (lohnston et al 



1992 ) , but is in reasonable agreement with th e inter feromet 



ric measurements of Frail, Goss & Whiteo ak ( 1994 ) and the 
new timing position of Wang et al. Our position for 
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PSR B 1754-24 has significantly smaller uncertainties than 
the catalogued position. The error in the latter is certainly 
much greate r than the beamwidth in the observations of 
Kijak et al. (1998), and can account for their 4.9 GHz non- 



detection of this flat-spectrum pulsar. The majority of the 
remaining positions are consistent, to within 3a, with posi- 
tions determined from pulsar timing. 

For the majority of sources in our sample, the off-pulse 
image showed no emission, extended or point-like, at or near 
the position of the pulsar. Sources of note are discussed in- 
dividually below: 

B0114+58: On-pulse data were corrupted for this source by 
hardware problems during the observations, and so no de- 
tection of the pulsar was made. The off-pulse image was un- 
corrupted, and shows no emission at the catalogued pulsar 
position. 

B0136+57: An unresolved source of flux density 2.6 ± 
0.3 mjy, approximately half of the observed pulsar flux den- 
sity, is seen at the pulsar's position in an off-pulse image. 
This source is ~ 90% linearly polarized, which is similar 
to the degree of polarisation seen for the pulsar (Gould & 
Lyne 1! 98). Thus the source probably corresponds to a com- 



ponent of the pulse-profile which was not properly gated out 
when the gate was set during observations. 
J0631+1036: Images of the region were badly corrupted by 
sidelobes from the source 4C+10.20 (flux density 2.5 Jy), 19' 
distant. The pulsar was not detected, its tabulated 1.4 GHz 
flux density of 0.8 mjy being below the sensitivity of the 
data. No other emission at or near the pulsar's position was 
detected, down to the sensitivity limit. 

B0656+14: No off-pulse emission is apparent at the pul- 
sar's position, but an extended, polarized sou rce is seen ~ 2' 
south of the pulsar, which Cordova et al. (1989) argue is 
possibly associated with the pulsar. We note that we see 
no connecting structure between the pulsar and this source, 
despite being more sensitive to extended emission than Cor- 
dova et al. (1989). It thus seems likely that this source is 



unrelated to the pulsar. We note that an X-ray PWN as- 
sociated with this source was claimed by Kawai & Tamura 
(1996 ), but has be en discredited by higher resolution data 
( Pecker et al. 199S| ). 

B 1356-60: The region around this pulsar is shown in Fig [j] 
- the pulsar lies on the western rim of a shell of radio emis- 
sion. This shell, which we designate G311.28+1.09, is ap- 
proximately circular, with a diameter ~9 arcmin and a flux 
density 0.04 ± 0.01 Jy at 1.4 GHz. In linear polarization 
there is significant confusion from other sources in the re- 
gion, and so it is not possible to determine whether the shell 
is polarized. Based on morphology alone, we thus consider 
G311.2 8+1.09 to be a possible new SNR. although further 



observa tions will be required to confirm this. However, given 
the 300- kyr characteristic age of the pulsar, it is unlikely that 
there is any physical association between PSR B1356-60 and 
G311.28+1.09. In an ungated image, no PWN is apparent 
around the pulsar, although the sensitivity to such a source 
is poor due to the presence of G311.28+1.09. 
B 1508-57: This pulsar is in a confused region of the Galactic 
Plane. We were forced to discard short u — v spacings in 
order to image the region, limiting the largest spatial scale 
to which these observations were sensitive to only 1'8. No 
off-pulse emission was observed, subject to this constraint. 
B1634~45: The pulse profile for this source reveals an inter- 



pulse separated in phase by 180 degrees from the main pulse, 
a result which has been confirmed by recent timing obser- 
vations (F. Crawford, private communication). Gating out 
both the pulse and interpulse reveals an unresolved off-pulse 
source at the pulsar's position of flux density 0.8 ± 0.4 mjy. 
This source is ~80% linearly polarized, similar to that mea- 
sured for the main pulsed component. Since no PWN has 
ever been observed to be so highly polarized, we think it un- 
likely that this off-pulse source corresponds to an extended 
nebula; it is more likely that this emission comes from the 
pulsar itself. While we cannot rule out an error in the gating 
hardware or software, we note that no other gated ATCA 
observations have shown such an effect (see SGJ99). Alter- 
natives are that there is a low-level bridge of pulsed emission 
connecting the two main components of the pulse profile 
(cf. PSR B1259-63; Manchester & Johnston 1995), or that 
the pulse profile contains an underlying unpulsed compo- 
nent (cf. PSR J0218+4232; Navarro et al. 1995). We are 
planning further ATCA observations of this source in order 
to distinguish between these possibilities. 
B1706-16: An unresolved off-pulse source of flux density 
4.0 ± 0.3 mjy is seen at the pulsar's position. This source is 
less than 15% linearly or circularl y polarized, but so is the 
pulsar itself (Gould & Lyne 1998). While the pulse profile 
is quite narrow and shows no evidence for an interpulse, 
the fact that the VLA gating is set on-line means that, as 
for PSR B0136+57, we are unable to rule out a component 
of the pulse profile as the source of this detection. As for 
PSR B1634-45 above, we plan to re-observe this pulsar with 
the ATCA in order to clarify this situation. 
B 17 18-35: This pulsar is in a complicated region, and suffers 
significant confusion from the nearby star-forming region 
NGC 6334 (e.g. Brooks & Whiteoak 2000). Gating shows 
the pulsar to be located at the center of a 4' radio nebula, 
G351. 70 +0.66, which is al so clearly visible in data from the 
MGPS (| reen et al. 1999] ) . This region has a distinct coun- 
terpart at 60 /im in IRAS data, and is probably thermal. 
B 1727-33: Two ultra-compa ct Hn regions are in the field, 
one of which, G354.19-0.06 ( pecker et al. 1994|) , has a flux 
density of 0.3 Jy and is located ~ 10' from the pulsar po- 
sition. The pulsar was not detected, its catalogued flux of 
2.9 mjy corresponding to a signal-to-noise of only 2a. 
B 1730-37: The sensitivity of the observations was reduced 
by the presence of PMN J1733-3722 (flux density 0.6 Jy), 
just 2' away. 

B1754S4: The la uncertainty in the right ascension of this 
pulsar was previously 14'; as discussed above, we have now 
greatly improved on this position. This more precise position 
puts the pulsar along the sa me line of sight as the large 
diffu se Hn region G5. 33+0. 08 (Lockman, Pisano & Howard 



199f). Emission from the latter is clearly seen in NVS S dat a 



and in the 90 cm image of Frail, Kassim & Weiler (1994) 
but is largely resolved out by our observations. 



4 DISCUSSION 

As discussed in §0, the parameters of the current survey were 
chosen to improve on the sensitivity of previous surveys, in 
particular that carried out by FS97. These factors are sum- 
marised in Fig ^, where the sensitivity of FS97's search is 
compared to that presented here. Only for PWN with radii 
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Figure 1. 1.4 GHz ATCA image of the region surrounding PSR B1356-60. Contours are at levels of 2, 4, 6, 8, 10 and 
12 mjy beam -1 , at a resolution of 46" X 26" (FWHM shown at lower right). The pulsar (which has not been gated out in 
this image) is marked with a cross, and can be seen sitting on the western rim of the shell G311.28+1.09. 



between (f.'4 and ~ 1" is FS97's search more sensitive than 
ours. Between 1" and 10", the current results are consider- 
ably (up to 100 times) more sensitive, while at scales smaller 
than Cf.'A and larger than 10", our data probe a parameter 
space to which FS97 were not sensitive at all. 

The results reported in showed that most of the 
sources in our sample had no detectable PWN associated 
with them. The exceptions were PSRs B0136+57, B1634-45 
and B1706-16, for which unresolved off-pulse sources were 
detected at the pulsar position. From the current data, we 
are unable to conclusively determine whether this emission 
corresponds to emission from the pulsar or from a compact 
PWN. While we plan to investigate these sources further, 
for the purposes of the present discussion we assume these 
observations to be non-detections, but with a sensitivity cor- 
responding to the flux density of the off-pulse source (rather 
than to the noise in the surrounding area of the image) . 

To quantify the significance of our non-detections, we 
follow FS97 in characterising a radio PWN's integrated lu- 
minosity, Lr, by 



Lb = eE er 



g s 



(1) 



where E erg s _1 is the associated pulsar's spin-down lu- 
minosity, and e is the fraction of E which goes into radio 
emission. 

Assuming a typical PWN spectral index a — —0.3 and 
integrating from 10 MHz to 100 GHz, the corresponding 
1.4 GHz flux density is 



Si.. 



2.1 x 10" 



eEs, 



mjy, 



(2) 



kpc is the distance to the pulsar and E = 
1 . We generally use distances from the pulsar 



where d 
10 34 £ 34 erg s 

catalogue ( Taylor et al. 1995), derived either from a pul- 
sar's dispersion measure flTaylor fc Cordes 1993; ) or from its 
kinematic distance based on H I measurements (e.g. Frail & 
Weisberg 1990). 

If we do not detect a PWN in our observations, we can 
potentially put an upper limit on Si. 4, and hence on La and 
e. However, as demonstrated in Fig ^, these limits depend on 
the angular size we expect for the PWN. As previously dis- 
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Radius of PWN (arcsec) 

Figure 2. Comparison of 5cr sensitivities of the 8.4 GHz PWN search of FS97 with the 1.4 GHz observations reported here. 
The observations of FS97 were not gated, and so generally are not sensitive to PWN on scales smaller than the resolution 
limit. For our VLA (ATCA) observations, we have adopted a typical spatial resolution of 12" (12"), a maximum spatial 
scale of 15' (4.'5), and a la sensitivity of 0.5 (0.15) mjy beam -1 . A spectral index a = —0.3 (S v oc u a ) has been assumed 
in converting 8.4 GHz results to 1.4 GHz. 



cussed by FS97 and SGJ99, PWN can in general be divided 
into two distinct classes (assuming that the pulsar is not in- 
side a SNR, which appears to be the case for all the sources 
in our sample) : those which are confined by the gas pressure 
of the ambient ISM ("static PWN"), and those confined by 
ram-pressure resulting from motion of the pulsar through 
the ISM ("bow-shock PWN"). 

Let us first consider the case of a static PWN. The 
bubble in the ISM driven by the pulsar wind will expand 
supersonic ally into the ambient medium, producing a PWN 
of radius ( Arons 1983j ) 



Rstm 



0.14 x 



/ E'j,jij 

V n 



1/5 



pc. 



(3) 



where n cm is the density of the ambient medium (as- 
sumed to be pure hydrogen), and t 3 kyr is the period for 
which the pulsar has been interacting with the ISM. In fur- 
ther discussion we assume that this age is given by the pul- 
sar's characteristic age, r c = P/2P. 

However a PWN is only static while i? s tatic > Vpsr, 
where Vpsr km s _1 is the pulsar's space velocity. Re- 
arranging Equation (3) of FS97 (and correcting for a missing 
factor of 4-7T in their results), we find that a PWN will be 
static for velocities and ambient densities for which: 



nVp'sn <> 4 x TO 9 E34/4 



(4) 



When this condition is not met, the pulsar has "overtaken" 
its own static PWN, and a bow-shock PWN results. The 
resulting PWN is much smaller than the static PWN, and 
has a size determined by a balance between the pressure of 
the relativistic pulsar wind and the ram-pressure resulting 
from the pulsar's motion, 



Rbo 



0.63 



E34 

nV PSR 



1/2 



pc, 



(•») 



where we have assumed that all the spin-down luminosity 
of the pulsar goes into the wind, and that the radius of a 
bow-shock PWN is 1.5 time s the radius at which the r am 
and wind pressures balance ( van Buren fc McCray 198£ ) . 

Thus for a given n and Vpsr, we can use Equation Q 
to determine whether a PWN has a static or a bow-shock 
morphology, and from this use either Equation (jij) or (|s|) to 
determine its radius, and hence its angular extent, #pwn- 

For a pulsar- gated observation with a 1.4 GHz RMS 
sensitivity of a mjy beam" 1 at a resolution 9'^ x 6g , sup- 
pose no off-pulse source is detected at the 5a level. There 
are three possible reasons for this non-detection: the PWN 
is unresolved and below the point-source sensitivity of the 
observations, the PWN is resolved and below the surface 



© 2000 RAS, MNRAS 000, [l[|] 



Limits on radio emission from pulsar wind nebulae 7 



brightness sensitivity limit, or the PWN is larger than the 
maximum spatial scale to which we are sensitive. We now 
consider what limits on e can be derived for each situation. 

If flpwN is smaller than the resolution limit, then Si. 4 < 
5cr and from Equation (^) the corresponding upper limit on 



e < e = 2.4 x 1CT 5 a (E 3 i/d 2 



(6) 



Note that if an off-pulse point source of 1.4 GHz flux density 
S p t is detected at the position of the pulsar (as was the case 
for three of the pulsars in our sample), then 5a should be 
replaced by S p t in this expression. 

If the PWN is extended, but smaller than the largest 
spatial scale to which our interferometer is sensitive, then 



e < 



(7) 



where we have assumed a Gaussian profile for the PWN. Fi- 
nally, if the PWN is larger than our observations can detect, 
we can put no limit on e. 

FS97 have argued, on the basis of their non-detections, 
that the lack of observable radio PWN around most pul- 
sars implies e ^ 10 -6 , two orders of magnitude less than 
for pulsars with detected radio PWN. However, in making 
this calculation they assumed that the ambient density was 
n — 1 cm -3 . Through Equation (Q), when combined with a 
reasonable velocity, this implied that all the pulsars they ob- 
served were powering bow-shock PWN, and the consequent 
high ram pressure ensured that these sources would largely 
be unresolved by their 0'.'8 resolution. 

However their assumed ambient density is not represen- 
tative of our current best understanding of the ISM. While 
the relative filling fractions are uncertain, available evidence 
supports a multi-phase ISM of which 90% by volume is a 
combination of a warm medium of density n — 0.3 cm~ 3 
and a hot ionised component of density n = 0.003 cm -3 
(see Ferriere 1998a for a recent discussion and overview). In 
the hot low-density component, we still expect most PWN 
to be bow shocks, but the ram pressure is greatly reduced 
and PWN will consequently be much more extended. Thus 
many PWN in this low density medium will be larger than 
the resolution limit of FS97, and through Equation ([?]), the 
limit on e much less stringent than claimed. 

In order to better consider the limits on e, we there- 
fore carry out the following calculation for 27 of the pulsars 
in Table |l| (PSR B 1706-44 is excluded as it has a candi- 
date PWN), and additionally for the 4 pulsars discussed by 
SGJ99. For each pulsar in our sample, we adopt possible 
densities of n = 0.3 cm -3 and n — 0.003 cm -3 . Ten of the 
pulsars in our sample have measured proper motions, and 
we use the correspo nding 3D space velocities determined by 
Cordes & Chernoff ( |l99g| ); one other pulsar (PSR B1055-52) 



has ha d a scintillation veloc ity determined for it (Johnston 



Nicastrs & Koribalski 1998). For the remaining sources we 
set Vpsr = 380 km s~ , corresponding to the m ean p ulsar 
velocity of the distribution of Cordes & Chernoff (1998). Us- 



ing Equation (^) we determine whether each corresponding 
PWN is bow-shock or static, and then consequently deter- 
mine fl pwN- Upper limits on e are then determined from 



Equati ons ( 6 ) and ( 7 ), and the results given in Table 2 . For 
the 16 sources in our sample which were also observed by 
FS97, we can make a similar calculation based on their re- 



sults (converting their data to 1.4 GHz assuming a spectral 
index a = —0.3); these revised values of e ma x are also listed 
in Table ^. Note that many PWN are unresolved for either 
choice of ambient density, and so have the same value of 
tmax in both cases. 

We can compare our upper limits to values of e for 
known PWN. In Table 3 of FS97, e is listed for the six pulsars 
then known to have detected radio PWN, to which we add 
PWN which have since been associated with PSR J0537- 
6910 (e = 5 x 10~ 4 ; Lazendic & Dickel 1998), PSR B0906-49 
(e = 2 x 10 -6 ; GSFJ98) and PSR J1811-1926 (e < 2 x 10 -3 ; 
Morsi & Reich 1987; Torii et al. 1999). For these nine 
sources, six measurements lie in the reasonably narrow range 
(1 - 5) x 10" 4 . Those sources lying outside this range are 
PSR B0833-45, for which it is unclear just what part of the 
surrounding SNR is pulsar-powered, PSR B0906-49, which 
GSFJ98 argue is substantially different from other PWN, 
and PSR J1811-1926, where the radio PWN is faint and 
has poorly constrained properties. Thus the data available 
suggest that a "typical" detectable PWN has e « 10~ 4 . 

FS97 argued that typical non-detections corresponded 
to e ma x ~ 2 x 10 -6 , significantly less than for detected PWN. 
However, it can be seen from the results in Table ^ that for 
n = 0.003 cm~ 3 some PWN become too large to be detected 
by their data, while for the remaining pulsars the limit rises 
to e max ~ 3 x 10~ 4 . Thus we argue that the observations 
of FS97 could have missed "normal" PWN around most of 
their sample if most of these sources are in low density re- 
gions, and hence their observations do not constrain pulsars 
which lack detectable PWN to be any different in their wind 
properties from those pulsars with observed PWN. 

On the other hand, the current observations can poten- 
tially provide constraining limits on e. We first consider the 
six young and energetic pulsars in our sample (i.e. the first 
six pulsars in Table ^|). These pulsars are defined approxi- 
mately by S34 J; 50 and tz Si 50, properties similar to those 
pulsars with detectable radio PWN. For either assumed ISM 
density, our data constrain these pulsars to have upper lim- 
its on e in the rang e (0.002 - 0.2) x 10~ 4 , significantly less 
than for pulsars with observed PWN. While this seems to 
imply genuinely low values of e, for n = 0.003 cm -3 , con- 
ditions are such that these pulsars have only recently over- 
taken their static nebulae. Lowering their velocities slightly, 
arguing that their actual ages are less than their character- 
istic ages, or accepting that realistically, the transition from 
static to bow-shock PWN does not happen instantaneously, 
it seems likely that these pulsars are still producing static 
nebulae, whose extents are much larger than for bow-shocks. 
In this case, the corresponding limits become e max ^> 10~ 4 , 
and are not constraining. 

We thus argue that our non-detections of PWN can be 
explained even if all young and energetic pulsars have sim- 
ilar wind properties. The difference between detectable and 
non-detectable PWN seems to be that detected PWN are 
either in dense regions of the ISM or in SNRs, in which 
there is sufficient external pressure to confine the pulsar 
wind and produce an observable PWN. However, pulsars 
with no PWN are in the low density phase of the ISM and 



so pr oduce unobservable "ghost remnants" ( Blandford et al 



1973; ). With the exception of PSR B0906-49 (GSFJ98), pul- 
sars with observed PWN also have associated SNRs, while 
all those young pulsars without PWN also have no associ- 
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Table 2. Upper limits on e = Lr/ E for PWN non-detections, including four non-detections from SGJ99. Pulsars are sorted by decreasing 
E\ 9 is the predicted angular size of a PWN for a given density. For upper limits on e derived from the results of FS97, U indicates that 
the predicted size of the PWN is too large to have been detected by their observations, while ". . ." indicates that a particular pulsar was 
not part of their sample. 



Pulsar 


E34 


*3 


d 


^PSR 




n = 0.3 cm 


3 




n = 0.003 cm 


-3 




(erg s _1 ) 


(kyr) 


(kpc) 


(km s — ) 


6 (") 


l°gio e 
(this paper) 


max 

(FS97) 


6 (") 


loglO e 
(this paper) 


nax 

(FS97 


B1823-13 


280 


21 


4.1 




5.01 


-6.3 


-5.4 


50.1 


-5.2 


U 


Jl 105-6107 


250 


63 


7.1 




2.76 


-5.8 




27.6 


-5.6 




B 1046-58 


200 


20 


3.0 




5.86 


-6.8 




58.6 


-5.3 




B 1610-50 


160 


7 


7.2 




2.16 


-5.8 




21.6 


-5.3 




B1727-33 


120 


26 


4.2 




3.19 


-5.2 


-5.3 


31.9 


-4.8 


U 


B1930+22 


75 


40 


12.1 




0.88 


-4.6 


-5.4 


8.8 


-4.6 


-3.4 


B0114+58 


22 


275 


2.1 




2.70 


-5.6 


-5.4 


27.0 


-5.1 


U 


J1835-1106 


18 


127 


3.1 




1.69 


-5.1 




16.9 


-5.0 




J0631+1036 


17 


44 


6.6 




0.78 


-4.0 




7.8 


-4.0 




B0740-28 


14 


157 


1.9 


276 


3.36 


-5.5 


-4.6 


33.6 


-4.9 


U 


B1508-57 


13 


298 


12.7 




0.35 


-4.2 




3.5 


-4.2 




B1356-60 


12 


319 


5.9 




0.72 


-4.7 




7.2 


-4.7 




B1634-45 


7.5 


590 


3.8 




0.88 


-5.1 


-5.3 


8.8 


-5.1 


-3.3 


B0611+22 


6.3 


89 


4.7 


212 


1.18 


-4.5 


-4.9 


11.8 


-4.5 


-2.9 


J0538+2817 


4.9 


619 


1.8 




1.54 


-5.2 




15.4 


-5.0 




B1718-35 


4.5 


176 


6.4 




0.41 


-3.5 




4.1 


-3.5 




B0355+54 


4.5 


563 


2.1 


210 


2.29 


-5.2 


-4.8 


22.9 


-4.7 


U 


B0540+23 


4.1 


253 


3.5 


348 


0.77 


-4.5 


-5.1 


7.7 


-4.5 


-3.3 


B 1754-24 


4.0 


285 


3.5 




0.71 


-5.0 




7.1 


-5.0 




B0656+14 


3.8 


111 


0.8 


331 


3.63 


-6.0 


-5.2 


36.3 


-5.1 


U 


B1719-37 


3.3 


345 


2.5 




0.89 


-4.6 


-5.4 


8.9 


-4.6 


-3.4 


B1821-19 


3.0 


574 


5.2 




0.41 


-4.4 


-4.7 


i.l 


-4.4 


-3.3 


B1055-52 


3.0 


535 


1.5 


440 


1.21 


-5.8 




12.1 


-5.7 




B2011+38 


2.9 


412 


13.1 




0.16 


-2.9 


-4.0 


1.6 


-2.9 


-3.4 


B0136+57 


2.1 


403 


2.9 


340 


0.69 


-4.3 


-4.6 


6.9 


-4.3 


-3.3 


B1449-64 


1.9 


1035 


1.8 


337 


1.04 


-5.5 




10.4 


-5.5 




B 1730-37 


1.5 


355 


3.5 




0.44 


-3.6 




4.4 


-3.6 




B1933+16 


0.51 


947 


7.9 


996 


0.04 


-3.0 




0.4 


-3.0 




B1706-16 


0.12 


1655 


1.3 


186 


0.59 


-3.5 




5.9 


-3.5 




B0736-40 


0.089 


3805 


2.1 


377 


0.21 


-3.8 




2.1 


-3.8 




B2148+63 


0.012 


36640 


13.6 




0.01 


-0.7 


-1.6 


0.1 


-0.7 


-1.6 



ated SNRs. We indeed expect SNRs to be undetectable in 
the hot component of the ISM (Kafatos et al. 1980; Gaensler 



Johniton 1995), consistent with our conclusion above that 



it is a low ambient density which causes a PWN around a 
young pulsar to be undetectable. 

A notable exception is PSR B1757-24, which is asso- 
ciated with both a SNR and a radio PWN despite having 



0.003 cm~ A (Frail fc Kulkarni 1991; Manchester et al 



199l|p|n this case, the pulsar is inferred to have a transvers e 
space velocity ~1500 km s _1 (Frail, Kassim & Weiler 1994); 



this not only supplies the necessary ram pressure to produce 
an observable PWN, but has caused the pulsar to overtake 
the shell of the associated SNR, re-energising the remnant 
with its passage. If it were not for the extreme velocity of the 
pulsar, neither the PWN nor the SNR would be detectable, 
as expected in a low density region. 

The relative numbers of young pulsars with detected 
and undetected SNRs/PWN suggests an approximate fill- 
ing fraction ~50% for the low density component of the 
ISM. This is somewhat more than recent estimates of 15- 



20% (Ferriere 1998b), but can be explained by the fact that 
we expect young pulsars to be preferentially located in low 
density regions produced by the powerful winds of their pro- 
genitors and by previous supernovae in the region. 

The majority of pulsars in our sample are considerably 



less energetic (£34 < 50) and older (tz <; 100) than pul- 
sars around which radio PWN have been observed. Values 
of e ma x for these pulsars are plotted in Fig ^, from which 
it can be seen that for either choice of ambient density, the 
distribution peaks around e m ax ~ 10 -5 . These pulsars have 
all long since overtaken their static PWN, and will have 
bow shock PWN for any sensible choice of n and Vpsr. The 
resulting size of such a PWN is not a strong function of 
n or Vpsr; to produce values of e ma x > 10 -4 , consistent 
with detected PWN, requires uniformly low space velocities, 
Vpsr 150 km s _1 , for these pulsars. However, only ~5% of 
pulsars are thought to be traveling at less than 150 km s _1 
( |Cordes fc Chcrnoff 1998] ). While the ages and distances we 
have used for these pulsars have their associated uncertain- 
ties, for bow-shock PWN values of e ma x are independent of 
age, and distances would have to uniformly be increased by 
a factor of three to shift the peak in e ma x to a value in agree- 
ment with that seen for detected PWN. The lower value of 
e we have derived for older pulsars is thus a result quite 
robust to the assumptions and uncertainties involved in its 
derivation, and we therefore argue that these pulsars have 
winds which are genuinely at least an order of magnitude 
less efficient at producing radio emission than the winds of 
young and energetic pulsars. 
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FS97 consider various reasons why older pulsars might 
appear to have lower values of e. Possibilities include: 

(i) that the PWN are resolved out by the observations; 

(ii) that an increasing fraction of E goes into pulsed X- 
rays and 7-rays (see Thompson et al. 1994); 

(iii) that their winds are dominated by Poynting flux 
rather than relativistic particles; 

(iv) that the injection spectrum of particles in the pulsar 
wind has shifted to higher energies. 

Our observations can conclusively rule out alternative 
(i), as we can detect PWN produced for almost all feasible 
values of n and Vpsr. While we cannot distinguish between 
the remaining three possibilities, we note that of detected 
PWN, that associated with the oldest pulsar, PSR B0906- 
49, also has the lowest value of e and the steepest spectral 
index (GSFJ98). Since the spectral index of a PW N is di- 
rectly related t o the spectrum of injected particles (Pacini 
&i Salvati 1973), this result tentatively suggests that the ef- 
ficiency of the wind in producing radio emission is related 
to the injection spectrum, and that alternative (iv) might 
then best explain the observations. 



CONCLUSIONS 

We have used pulsar-gating at 1.4 GHz to search for radio 
PWN around 27 pulsars. Our search was up to 100 times 
more sensitive than the only other comparable survey, and 
was carried out on spatial scales corresponding to a much 



wider range of ambient densities and pulsar velocities. In- 
cluding data from previous work by SGJ99, non-detections 
towards 28 pulsars, plus inconclusive results in three other 
cases, have allowed us to determine upper limits on e, the 
fraction of a pulsar's spin-down luminosity which goes into 
producing radio emission from a PWN. 

We find that the data are consistent with virtually all 
young energetic pulsars having e ~ 10 -4 . The lack of PWN 
around ~50% of young pulsars can be explained if they 
are in low ambient densities (0.003 cm -3 ), consistent with 
the absence of associated supernova remnants around these 
sources. 

For older pulsars, any reasonable choice of ambient den- 
sity and pulsar velocity results in upper limits on the wind 
efficiency e < 10 -5 , ten times less than for young pulsars. 
Thus pulsars seem to become less efficient at producing ra- 
dio wind nebulae as they age; we speculate that this result is 
due to the spectrum of injected relativistic particles steep- 
ening in older pulsars. This possibility can be tested through 
X-ray observations towards such pulsars - it is likely that 
Chandra will make many new detections of X-ray PWN, and 
through consequent imaging spectroscopy, we may finally be 
able to probe the winds around older pulsars. 

Of the ~55 non-recycled pulsars with E > 3 x 
10 34 erg s _1 , almost all have now been searched for associ- 
ated radio PWN down to a good sensitivity. If those sources 
with no detectable PWN are indeed in low density regions 
of the ISM, it seems unlikely that we will ever find radio 
PWN around them with current telescopes. For example, if 
PSR B1046-58 is powering a static PWN with e = 10" 4 , 
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the resulting radio nebula would be 20' across with a flux 
density at 1.4 GHz of 0.5 Jy. To detect this source would 
require a = 0.3 mjy arcmin" 2 , which is generally below the 
confusion limit for instruments capable of imaging sources 
this large. We might have to wait for the large increase in 
sensitivity promised by the Square Kilometre Array in order 
to make further progress. 
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